Electron-phonon effects on the Raman spectrum in MgB 2 
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The anomalous features of the Raman spectroscopy measurement in MgB2 represent a still un- 
resolved puzzle. In particular highly debated are the origin of the huge E2 g phonon linewidth, the 
nature of the low energy (oj < u>E 2g ) background and the evolution of the Raman spectra with Al 
doping. In this paper we compute the self-energy of the i?2 9 phonon mode in a fully self-consistent 
way taking into account electron-phonon effects on the electronic properties. We show that all the 
anomalous features can be naturally understood in a framework where the whole electron-phonon 
spectrum a 2 F(w) gives rise to significant damping processes for the electronic excitations and con- 
sequently for the E2 g phonon itself. The two-peak structure as function of the Al doping is ascribed 
to finite bandwidth effects arising as the Fermi level approaches the a band edge. 

PACS numbers: 74.70.Ad, 74.25.Kc, 63.20.Kr, 78.30.-j 



There is nowadays a general consensus that the 
electron-phonon (el-ph) interaction is responsible for the 
superconductivity in MgB 2 El II II II II HI Inter- 
esting peculiarities make however this compound quite 
unique and different from the conventional low-T c super- 
conductors. Most notable are the two-bands/two-gaps 
phenomenology 0, [1 H H El EI j the fact that a large 
part of the electron-phonon coupling is concentrated in 
just one phonon mode Ei g 0, H U the possibility of 
spanning with Al doping a wide range of hole doping of 
the a bands E2 ■ Although several features are now well 
understood within the context of the two-gaps scenario, 
the experimental overview of the Raman spectroscopy of 
this compound stil l rep resents a highly puzzling anomaly 
[3 El El 13 El EI El The main open questions on 
this point regard in particular the origin of the anoma- 
lously large linewidth of the E^g phonon mode, the na- 
ture of the background Raman signal at low frequencies, 
and the evolution of the Raman spectra upon Al doping, 
which shows a transfer of spectral weight between differ- 
ent peaks more than a continuous hardening of the Ei g 
phonon mode. 

In this Letter we show that all these unconventional 
features can be explained within the electron-phonon 
scenario by the coexistence of a strongly coupled (E2 g ) 
phonon mode at q = and of a wide background 
electron-phonon scattering with other phonon branches 
and different momenta. This gives rise to a Fano-like 
phenomenology where the shape and linewidth of the 
Ei g mode is mainly related to the electron-phonon decay 
processes triggered by the other modes and by impurity 
scattering. Similar effects give rise to a low energy peak, 
equivalent to the Drude-like feature in the optical con- 
ductivity, which is responsible for the low energy signal. 
We show also that the transfer of spectral weight between 
different Raman structures as function of the Al amount 
can be a direct consequence of the vanishing of the Fermi 
energy. 



The origin of the broad phonon linewidth in Ra- 
man spectroscopy has been widely discussed in liter- 
ature |il El El IH E3 Common feeling on this 
subject is that the linewidth of a phonon mode 7 q re- 
flects the strength of the electron-phonon coupling A q of 
this specific mode, in the spirit of the Allen's formula 
7 q = 27riV(0) A q w q 23], where iV(0) is the electron den- 
sity of states (DOS) coupled with the phonon mode and 
cij q is the phonon frequency. First-principle calculations 
however indicate that the electron-phonon coupling alone 
of the Ei g phonon mode is not sufficient to explain the 
large experimental Raman linewidth |22j. In addition, 
the direct employment of the Allen's formula is strongly 
questioned in the q = case, relevant for the Raman 
spectroscopy, where a more careful analysis shows that 
the damping of a q = phonon mode scattering with a 
non-interacting electronic system is strictly zero • 

In contrast to the above scenario, a Quantum Field 
Theory analysis specifically addressed to the q = case 
has been developed in Refs. [H H3- Although these 
studies were mainly aimed to investigate phonon anoma- 
lies in the superconducting state, this approach can be 
as well employed in the normal state. Along this line, 
for instance, Marsiglio et al. remarked [H that the 
imaginary part of a q = phonon self-energy in the weak- 
coupling limit is simply 
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where the impurity scattering rate ri mp was assumed to 
be the only damping source of the electronic propagator: 
G(k, lo) = l/[w - e k + iT imp ]. 

Eq. explicitly shows that the damping processes of 
a q = phonon are triggered by corresponding damp- 
ing processes of the electronic charge response which 
screens the phonon, whereas the el-ph coupling A q= o of 
this particular mode rules the magnitude of these effects. 
This results holds true also when the main source of the 
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the electronic excitation damping is the electron-phonon 
interaction itself. In this case it is also important to 
note that in principle all the phonon modes, and not 
only the specific A q= o mode, contribute to the electronic 
damping. Indeed, as we are going to discuss, the imag- 
inary part II" (u) of the phonon self-energy is roughly 
related to the integral up to the energy u> of the imag- 
inary part of the electronic self-energy, 7(q = 0,lu) oc 
A q= o Jq du'TUj'). For optical phonons close to the top 
of the phonon spectrum we have f" duj'T(u') oc Aw, 
A being the total electron-phonon coupling constant, 
so that in intermediate-strongly coupled systems the 
phonon linewidth can be of the same order of the phonon 
frequency itself. In conventional low-T c superconduc- 
tors however the total electron-phonon coupling A ~ 1 
is spread over several phonon branches and q modes, 
A = X)qiy^q,^j so that the contribution from a single 
mode is quite small. Typical values of A q= o,i/ for the A g 
phonon modes in cuprates are for instance A q= o,i, < 0.034 
|26j , resulting in experimental Raman linewidths of few 
meV [53 ■ The opposite extreme case of only one strong 
coupled q = phonon mode loq carrying the whole 
electron-phonon coupling would also predict extremely 
small phonon linewidth since, although A q= o ~ 1, the 
amount of electron-phonon coupling smaller than loq is 
practically negligible and J Q ° duiT(nj) as 0. 



part reads: 
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where N c is the number of sampling point in the Bril- 
louin zone, A(k, cj) = (l/7r)Im[l/(w — £k — ^( w ))] is the 
electron spectral function and N s = 2 and Nb = 2 rep- 
resent respectively the spin degeneracy and the cr-band 
degeneracy. Note that Eq. J2J does not account the 
double degeneracy of the two E 2g modes which on the 
other hand contribute to the superconducting pairing. 
A similar expression is obtained for the real part of the 
phonon self-energy, also attainable from the Kramers- 
Kronig relations. The electron-phonon matrix elements 
9k,E 2g can be estimated from the E 2g deformation poten- 
tial D E2g ~ 12 ev/AHHE3,E3, 

which is in good ap- 
proximation k-independent close to the Fermi level, and 
from evaluating the zero point motion lattice displace- 
ment gE 2g — Dwm. y/ (u 2 ) = 0.39 eV, taking into account 
anharmonic effects [32j . 

We finally evaluate the phonon spectral function BUS) 
from the imaginary part of the phonon propagator 
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Along this scenario MgB2 presents very peculiar char- 
acteristics since it presents at the same time two (de- 
generate) strong coupled q = E 2g phonon modes and 
a relevant fraction of the total electron-phonon couglin 
strength spread over other different modes 0, 0, |T 
We employ a fully self-consistent many-body approach 
to investigate the normal state E 2g Raman spectrum 
of MgB2. In more explicit way we first solve itera- 
tively the Marsiglio-Schossmann-Carbotte equations 
with the Eliashberg spectral function o? a FUS) obtained 
by first-principle calculations [29| (Fig. QJi) to evaluate 
the real-axis electronic self-energy £(w) of the a bands. 
The Eliashberg function c%F{w) = a 2 aiJ FUj) + ofc^FUj) 
describes the total electron-phonon scattering of the a 
electrons with both the a ad 7r bands. Note that the 
a 2 F(u>) extracted from first-principle techniques repre- 
sents the el-ph spectral function where phonons are al- 
ready renormalized. In order to better compare with 
the experiments, we include also explicitly possible ef- 
fects of impurity disorder. We schematize the u-bands 
as two degenerate bands with constant DOS, namely: 
X; k G(k,w) = N a (0) J deG(e,uj), where N a (0) represents 
the cr-band electron density of states per spin and per 
band, N a (0) ~ 0.075 states/(eV • cell) 0,13]. The so- 
obtained electronic Green's function is then employed as 
input to calculate the full frequency dependence of the 
q = phonon self-energy HUj) |25j, whose imaginary 



where ^E 2g is the unrenormalized phonon frequency, 
while the renormalized phonon frequency WE 2g and the 
phonon linewidth jE 2g are simply related to the real and 
imaginary parts of n(w), 



2 

^e 2 



1E 2 



E 2g 



2Q E2g Il'(Lj E2g ), 



' 2 (^E 2g /uJE 2g )ll"(0JE 2 



(4) 
(5) 



We assume the unrenormalized Ei g phonon frequency 
^E 2g — 100 meV, close to the top of the phonon spec- 
trum. 

In Fig. \3p,c we plot the imaginary part of the phonon 
self-energy, and the fully renormalized phonon spectrum. 
The filled symbols mark the values of imaginary parts 
of the phonon self-energy evaluated at the renormalized 
phonon frequency u>e 2 obtained from the self-consistent 
solution of Eq. (0J. At low temperature and in the ab- 
sence of impurity scattering, the imaginary part of the 
phonon self-energy shows a monotonic behavior which 
reflects the corresponding increasing of the electronic 
damping processes (Fig. Ell). Already in this case we 
predict a relevant value of T\"UoE 2g ) and of the phonon 
linewidth jE 2g — 6 meV, which is essentially only due 
to the spectral weight of a^FUj) for w < uJE 2g - Things 
are even more drastic when finite temperature effects or 
impurity scattering are taken into account. As is well 
known the imaginary part of the electronic self-energy 
starts from a finite value at lu — which is reflected in a 
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FIG. 2: (a) Phonon spectral function for T = 
50, 100, 150, 200, 250, 300 K; (b) temperature dependence of 
the renormalized phonon frequency UE 2g ( c ) temperature de- 
pendence of the phonon linewidth 7_E 2g . We considered here 
a small amount of impurity concentration corresponding to 
Time = 5.5 meV. 



FIG. 1: (color online) (a) Electron-phonon spectral function 
o^,F{w) (from Ref. 29]). The dashed line represents the 
unrenormalized phonon frequency ^E 2g — 100 meV. Panels 
b,c,d show respectively the imaginary part of the phonon self- 
energy, the phonon spectral function and the imaginary part 
of the electronic self-energy for three representative cases. 



sudden increase of |II"(w)| and in a significant broaden- 
ing of the phonon spectrum. Note also the appearance of 
an incoherent background for ui < we, and of a broad 
shoulder at very low temperature, in striking agreement 
with the Raman experimental data. Similar effects were 
pointed out in Refs. [13, 13 in the context of the elec- 
tronic Raman scattering. We would like to stress that 
the sharp increase of |II"(w)| and the onset of the low 
temperature shoulder in the phonon spectrum are intrin- 
sic features of the charge response function. As a matter 
of fact, as noted in Ref. |25| . at a first approximation the 
ratio —U"(lo)/uj is qualitatively similar to the real part 
of the optical conductivity a' (ui) . The low energy sudden 
enhancement of |II"(w)| and the corresponding shoulder 
in the phonon spectrum, which arise from the finite value 
of E"(w = 0), have thus a strict connection with the ap- 
pearance of a Drude-like peak in a'(u>). Along this line, 
the experimental observation of a si gnifi cant Drude-like 
scattering rate 1/r ~ 9 — 37 meV |35l l36| and of the 
Raman shoulder at low temperatures T ~ 40 — 45 K 
[lij l points out the actual presence of a small amount of 
impurity scattering. From the comparison between our 
results and experimental data we estimate an impurity 



scattering rate r; mp ~ 5.5 meV. 

In Fig. |2t we show the evolution of the phonon spectral 
function B(lo) at function of the temperature. Increasing 
temperature leads thus not only to a smearing of the main 
phonon peak by increasing "fE 2g but also to a smearing 
of the low energy structure, in fair agreement with the 
experiments. In Fig. |2bj c we show also the temperature 
dependence of the renormalized phonon frequency oje 29 
and of the phonon linewidth "fE 2g ■ The agreement with 
the experimental data is once more remarkable for "fE 2g 
(lij. On the other hand our analysis underestimates the 
phonon frequency uje 2!I ^ 55 — 70 meV (to be compared 
with the u>E 2g ~ 78 meV from Raman spectroscopy) and 
predicts a phonon softening as a function of T which has 
not been observed. We remind however that anharmonic 
effects are not taken into account here, and that they 
would lead to a temperature dependent hardening of the 
E2 g phonon frequency. Since the anharmonic effects were 
shown to be related to the eletron-phonon interaction 
itself |30j ] , a correct evaluation of them would require once 
more the inclusion in a consistent way of the damping 
electronic processes. 

As a last point of our analysis we address the evolu- 
tion of the Raman spectra as function of the Al doping. 
Ab-initio calculations predict that the electron doping 
induced by the Al content would fill the two a bands, 
leading to a reduction of the electron-phonon coupling. 
This is expected to result in a steady hardening of the E 2g 
phonon frequency |3l| . Actual Raman measurements are 
however in substantially disagreement with this picture, 
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FIG. 3: Evolution of the phonon spectral function by de- 
creasing the a band Fermi energy for T — 40 K and Fi m p = 0. 
Inset: real part of the phonon self-energy (solid lines) for dif- 
ferent Ef, from top to bottom: Ef = 10, 20, 30, oo meV. The 
crossing with the dashed line determines the phonon peaks. 

showing that the phonon hardening occurs via a transfer 
of spectral weight from the Ei g peak to another, at the 
moment unknown, new high energy structure GHH. In 
the following we provide a qualitative scenario where this 
behavior is naturally explained in terms of finite band- 
width effects which arise when the Fermi level approaches 
the a band edge. 

A natural modelization to investigate this scenario 
would be to consider a semi-infinite band system 
J^^de — > where the Fermi energy Ef plays the 

role of an energy cut-off for the electronic excitations. 
However the correct treatment of the charge conserva- 
tion and of the asymptotic behavior at w — ► — oo in 
the real axis phonon self-energy present different nu- 
merical difficulties which make this approach unafford- 
able. For a qualitative insight we consider thus in the 
following a simpler model with a symmetric band /jfjj 
where these numerical problems are overcome. Since the 
main role is played by the presence of a finite cut-off Ep , 
this model is expected to shed a qualitative light also 
on more realistic cases. In Fig. [3] we show the evolu- 
tion of the phonon spectral function for Ef approaching 
zero for T — 40 K and Fi mp = 0. For Ep higher than 
60 meV (not shown here) the phonon spectrum presents 
only a one-peak structure corresponding to the renor- 
malized phonon frequency. As soon as 2Ep decreases, 
however, this structure rapidly looses spectral weight in 
favor of the unrenormalized phonon frequency £1 E 2g ac- 
companied by a narrowing of the phonon linewidth. This 
behavior can be understood by looking at the real part 
of the phonon self-energy: for Ef — > oo it has a negative 
minimum at roughly the maximum energy of a^F (ui), 



and then it vanishes very slowly for u> — > oo. Peaks 
in the phonon spectrum are qualitatively determined by 
Eq. (@J, corresponding to the crossing of W(u>) with the 
dashed line [(w 2 - f}| 2fl )/2Q B2 J in inset of Fig. [2 The 
finite band edge Ep is reflected in a sharp drop of II' (w) 
due to the cut-off in the particle-hole excitation, with 
a width roughly given by Ef. Decreasing Ef the de- 
termination of u>E 2g according Eq. has an abrupt 
jump from a renormalized (II' (oj) ^ 0) to the unrenor- 
malized (H'(uj) ~ 0) value. The crossover between these 
two regimes occurs roughly at -Ef — 20 meV where Eq. 
Q has three solutions, corresponding to the two peaks 
and the middle minimum in the phonon spectrum. 

We have been unable to show evidence of this behavior 
even at room temperature and in the presence of impurity 
scattering r; mp = 6 meV, where the large broadening of 
the phonon spectrum does not permit to resolve a two 
peak structure. One has to take into in account however 
that approaching the band edge also the electron DOS 
and the impurity scattering rate Fj mp oc N(0) vanishes. 
The resemblance of this framework with the experimental 
results [3, UJj is anyway striking suggesting that finite 
bandwidth effects are actually the natural explanation of 
the evolution of the Raman spectrum with Al doping. 

In conclusion in this paper we have investigated the 
phonon Raman spectroscopy data by computing the 
phonon self-energy of the E^ g mode. We show that the 
anomalous features of the Raman measurements, namely 
the huge phonon linewidth, the low energy background, 
the two-peak structure as function of the Al doping, can 
be naturally explained by the interplay of the E2 g phonon 
mode with the whole electron-phonon spectrum which 
gives rise to damping processes in the electronic excita- 
tion and in the Ei g mode itself. 
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